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Abstract 

We show that recent experiments based on Josephson junctions, 
SQUIDS, and coupled Josephson qubits have a cosmological inter- 
pretation in terms of axionic dark matter physics, in the sense that 
they allow for analogue simulation of early- universe axion physics. We 
discuss new experimental setups in which SQUID-like axionic inter- 
actions in a resonant Josephson junction environment can be tested, 
similar in nature to recent experiments that test for quantum entangle- 
ment of two coupled Josephson qubits. The parameter values relevant 
for early-universe axion cosmology are accessible with present day's 
achievements in nanotechnology. 
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1 Introduction 



What do advanced Josephson junction technologies, SQUIDs, coupled Joseph- 
son qubits and related superconducting devices used in nanotechnology have 
in common with the problem of dark matter in the early universe? A lot more 
than might seem obvious at first sight. One of the major candidates for dark 
matter in the universe is the axion. The equation of motion of the axion 
misalignment angle and that of the phase difference in a Josephson junc- 
tion are identical if the symbols in the mathematical equations are properly 
re-interpreted. This allows for analogue simulation of early-universe physics 
using superconducting electronic devices such as SQUIDS and Josephson 
junctions. It also allows for new experimental setups that test axionic in- 
teraction strengths in a Josephson junction environment, similar in nature 
to recent experiments that test for quantum entanglement of two coupled 
Josephson qubits. We show in this paper that the parameter values rel- 
evant for early-universe axion cosmology are accessible with present day's 
achievements in nanotechnology. Moreover, we will discuss novel types of 
SQUID-like interaction states by which axionic dark matter may couple into 
a given resonant Josephson junction environment. This may pave the way 
for novel dark matter detectors in the future. 

The research fields of superconductivity and cosmology are normally pro- 
ceeding independently of each other, and the two scientific communities don't 
know each other — dealing apparently with very different subject areas. But 
a look at the equations of motions of Josephson junctions, SQUIDS, coupled 
Josephson qubits and similar superconducting devices used in nanotechnol- 
ogy on the one hand and axionic dark matter on the other indicates that it 
makes sense to think about common approaches in both areas. The equa- 
tions of motion are very much the same (with a suitable re-interpretation 
of the currents involved) and hence it makes sense to develop common ap- 
proaches and to translate known results from one area (Josephson junctions) 
into possible results and phenomena for the other area (axion cosmology). 
This will be worked out in this paper. 

There are currently two major candidates for dark matter in the universe, 
with extensive experimental searches for both of them, WIMPS (weakly in- 
teracting massive particles) [TJ |2j [3] and axions [U El EJ [7J. In contrast to 
WIMPS, axions are very light particles. In spite of their small mass, they 
behave like a very cold quantum gas [5|. There are several experiments that 
search for axions in the laboratory, using e.g. cavities and strong magnetic 
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fields, which trigger the decay of axions into two microwave photons J5JEIIH]- 
These microwave photons are detectable if the cavity resonates with the ax- 
ion mass. Searches have been unsuccessful so far, but one needs to scan a 
huge range of cavity frequencies. For this purpose, SQUID amplifiers with 
very low noise levels are a very useful technological tool to reduce the noise 
level and to improve the scanning efficiency [8J. These applications of low- 
level noise SQUID amplifiers for axion searches are very important from an 
efficiency point of view for cavity experiments, but different from the fun- 
damental physics analogies between axions and Josephson junctions that we 
want to emphasize in this paper. Novel ideas of experimental axion detec- 
tion have recently also been presented in [HJ [TO] , based on cold molecules as 
suitable detectors [H] and analysis of X-rays from the sun [TU] . Quasi-axionic 
particles play also an important role in topological insulators, new materials 
with exotic metallic states on their surfaces [TTj [T2| [T3~| [H] . All this illus- 
trates that it does make sense to look at axion physics in a much broader 
context than within the original model, which was dealing with the strong 
CP problem in the standard model of elementary particle physics |15j . 

The axion is described by a phase angle, the axion misalignment angle. 
A Josephson junction is also characterized by a phase, namely the phase dif- 
ference of the macroscopic wave function describing the two superconducting 
electrodes of the junction. Our major motivation in this paper is the fact that 
the equation of motion of the axion misalignment angle is identical to that of 
the phase difference of a resistively shunted Josephson junction, with a suit- 
able re-interpretation of the symbols involved. As a first step, this opens up 
the theoretical possibility to connect both fields, and to make analogue ex- 
periments simulating axion cosmology using superconducting devices in the 
laboratory. Moreover, and more importantly, this novel approach also opens 
up the possibility to test for interaction strengths of incoming present-day 
axionic dark matter in a given resonant Josephson junction environment. 
The principal idea is that axions may briefly form SQUID-like interaction 
states in a resonant Josephson environment, before decaying into microwave 
photons. There are no cosmological constraints on these types of dark matter 
interactions, since almost all of the matter of the universe is not in the form 
of Josephson junctions. The only way to either confirm or refute such a hy- 
pothesis is by doing laboratory experiments. Suggested future experiments of 
this type can be easily performed in the laboratory and may ultimately open 
up the route for novel methods of dark matter detection based on modern 
nanotechnology. 
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This paper is organized as follows. In section 2 we point out the equiva- 
lence of the field equations of axions and Josephson junctions. In section 3 we 
show that there is quantitative agreement of the relevant parameters. In sec- 
tion 4 we investigate SQUID-like interaction states of axions and Josephson 
junctions. Finally, in section 5 the axionic Josephson effect is discussed. 

2 Comparing the equations of axion- and Joseph- 
son junction physics 

Let's compare the mathematics underlying both axions and Josephson junc- 
tions. Consider an axion field a = f a 9, where 9 is the axion misalignment 
angle and f a is the axion coupling constant. In the early universe, described 
by a Robertson Walker metric, the equation of motion of the axion misalign- 
ment angle 9 is 

2 4 
777 p 

9 + 3H9 + -?—sm9 = 0, (1) 
nr 

neglecting spatial gradients. Here H is the Hubble parameter and m a denotes 
the axion mass. The forcing term sin 9 is due to QCD instanton effects. In a 
mechanical analogue, the above equation is that of a pendulum in a constant 
gravitational field with some friction determined by H. 

When electric and magnetic fields E and B are present as well, the axion 
couples as follows: 

9 + 3H9 + ^ sin 9 = ?l±-c 3 e 2 EB. (2) 

g 1 is a model-dependent dimensionless coupling constant (g 1 = —0.97 for 
KSVZ axions, g 1 = 0.36 for DFSZ axions). The typical parameter ranges 
that are allowed for dark matter axions are [5j |6] 

6 • 10- 6 eV < m a c 2 < 2 • lO^el/ (3) 

and 

3 • 10 9 GeV <f a < 10 12 GeV. (4) 

The product m a c 2 f a is of the order m a c 2 f a ~ 6 • 10 15 (ey) 2 . 

Let's now look at resistively shunted Josephson junctions (RSJs) [T6| ITT]. 
In the 'tilted-washboard' model the phase difference S of the macroscopic 
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wave function of the two superconductors satisfies 



Here R is the shunt resistance, C is the capacity of the junction, and J c is 
the critical current of the junction. The frequency 



is the plasma frequency of the Josephson junction. The product 

Q := ujRC (7) 

is the quality factor of the junction. 

If a bias current I is applied to the junction by maintaining a voltage 
difference V between the two superconducting electrodes, then the equation 
of motion is 

■■ 1 • 2el c . . 2e T 

Remarkably, the equations of motion of axions and of RSJs are identical 
provided the following identifications are made in eqs. (j2J) and (jSJ): 

3H - w < 9 > 



m^c 4 2el r 



a 



h c 

2e 



(10) 



^EB = -J. (11) 

An interesting consequence of this observation is the fact that it is possible 
to make analogue experiments with RSJs that simulate axion cosmology 
in the laboratory. To simulate an axion in a certain era of cosmological 
evolution, one builds up an RS J with parameters R, C, I c , I fixed by eqs. (jHJ)- 
( jTTj) . The left-hand side are cosmological parameters, the right-hand side 
is nanotechnological engineering. Essentially the inverse Hubble parameter 
H^ 1 (the age of the universe) fixes the product RC, the axion mass fixes the 
critical current I c and the axion coupling to external electromagnetic fields 
EB/f% is represented by the strength of the bias current /. 
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3 Coincidence of axion parameters with those 
of superconducting devices 

We now show that there is not only qualitative, but also quantitative co- 
incidence between axion and Josephson junction physics. As worked out 
below, the numerical values of the parameters for typical axionic dark mat- 
ter physics and for typical Josephson junction experiments have similar or- 
der of magnitude. Let us start with the experiments performed by Koch, 
Van Harlingen and Clarke in [18J. They built up four different samples of 
Josephson junctions with parameters values in the range R ~ 0.075 — 0.77fi, 
C ~ 0.5 - 0.81pF, I e ~ 0.32 - 1.53mA. According to eqs. ©-(HI]), these ex- 
periments of Koch et al. simulate axion-like particles in a very early universe 
whose age is of the order H" 1 = 3RC ~ 10~ 13 — 10 _12 s. The axion mass is 
in the range 0.98 ■ 10~ 3 — 1.58 • 10~ 3 eV. This simulated axion mass is just at 
the upper end of what is expected for dark matter axions, see eq. ([3]). 

The recent experiment by Steffen et al. [19J, which tests for quantum 
entanglement of coupled Josephson qubits, is classically well described by an 
RSJ model where the product 3RC is of the order H~ l = 3RC = 1.2 • 10 _6 s 
[2D] . This experiment thus simulates the dynamics of weakly coupled axions 
after the QCD phase transition in the early universe, which took place during 
an epoch where H~ x ~ 10~ 8 s. It corresponds to an axion mass of m a c 2 = 
hu = 3.3 ■ 10 _5 eV, much smaller than for the Koch et al. experiment [18] 
but within the range expected for dark matter axions, see eq. ([3]). 

It is interesting to see that recent experiments [TDJ, (2DJ [211 [221 [231 121] 
simulate axions with masses in the entire range of what is interesting from 
a dark matter point of view. The experiments of Penttillae et al. [21] , deal- 
ing with superconductor-insulator phase transitions, correspond to m a c 2 = 
1.32 • 10~ 4 eV. Nagel et al. [22] deal with negative absolute resistance ef- 
fects in Josephson junctions, these experiments formally have m a c 2 = 2.83 ■ 
10 _5 eV. Superconducting atomic contacts [23] correspond to even smaller 
axion masses, namely m a = 6.7 • 10~ 6 eV, at the lower end of what is allowed 
in eq. ([3]). Two-dimensional Josephson arrays, as built up in [24J, correspond 
to arrays of coupled axions with m a c 2 in the range 6.62 ■ 10~ 5 — 1.52 • 10 _4 eV. 
All these recent experiments are within the range of axion masses that are of 
interest from a dark matter point of view. They can thus be regarded as re- 
alistic analogue experiments simulating an axionic dark matter environment 
in the early universe. 
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The main idea of this paper is to go a step further. If experiments of the 
type mentioned above simulate a kind of realistic axionic dynamics, could 
these experiments then be used to detect current day incoming dark matter 
axions, by triggering their decay in a resonant environment? This will be 
worked out in the following sections. 

4 Possible interaction mechanism between ax- 
ions and Josephson junctions 

Consider a RSJ which has a plasma frequency u p = J close to the dark- 
matter axion mass m a , according to eq. fflUj) . and that is driven by an external 
bias current / > I c . Free axion quanta correspond to small nearly-harmonic 
oscillations of the misalignment angle 8, in accordance with eq. ([I]). Suppose 
such an axion enters a Josephson junction that has similar parameters as 
the entering axion. Then this basically means we have a system of two 
Josephson junctions, the second one represented by the entering axion. The 
axion is very cold and its potential is given by a tilted washboard potential. 
Thus its phase variable may couple into a given Josephson environment in a 
SQUID-like configuration. This is illustrated in Fig. 1. 

For any Josephson junction, a given numerical value of a phase difference 
A(p on its own does not have physical meaning since it is not a gauge-invariant 
quantity (see, e.g., [IT]). Rather, what is of physical relevance is the gauge- 
invariant phase difference 

5 = Aw - — / Ads, (12) 
$o J 

where the integration path is from one electrode of the RSJ to the other (A is 
the vector potential under consideration and $o = ^ is the flux quantum). 
The term J Ads has profound consequences for the physics of Josephson 
junctions and SQUIDS if magnetic fields are applied. 

The gauge- invariant phase of the Josephson-axion-SQUID of Fig. 1 must 
be single-valued. Putting a closed integration path that passes the weak 
links as well as the interior of the superconductor, and using the same line 
of arguments as for ordinary SQUIDS [17J, one obtains 

5-6 = 2tt— mod 2vr. (13) 

<±> 
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Figure 1: A Josephson junction (junction 1) and an axion (formally rep- 
resented by junction 2) interacting in a SQUID -like configuration. The 
supercurrent / flowing through junction 1 simulates to the synchronized ax- 
ion the formal existence of a very large product E ■ B given by eq. f[TT|) . Thus 
the axion immediately decays into two microwave photons. 
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Here $ is the magnetic flux enclosed by the SQUID. If the flux $ is close 
to zero or given by integer multiples of $o> then the above condition simply 
implies 

5 = 9, (14) 

i.e. the two phases synchronize. 

Classically, the coupling between two Josephson junctions with shunt 
resistance R, capacity C, and inductivity L is described by the following 
coupled differential equations [20J: 

5 + -^-5 + u 2 sm5 = 1x {9-5) + —{5 + 2<kM 1 ) 

+ ^ + uj2sin0 = 7*(*-0) + 7^r(0 + 27rAf 2 ) 

(15) 

Here 5 is the phase difference of the first junction, 9 that of the second 
junction. Mj = $j/$o is the normalized flux enclosed by junction i (i = 1, 2), 
and 7 X = C x /C is a small dimensionless coupling constant, assuming both 
junctions are capacitively coupled by a capacity C x . For example, in the 
experiments of Steffen et al. dealing with coupled Josephson qubits [19] one 
has '-fx = 2.3 • 10~ 3 . Typically, the damping term proportional to 5 and 9 is 
neglected in the theoretical treatment of these types of experiments [20]. The 
above classical equations of motion describe quite well the experimentally 
observed phenomena [TjJl [20] • 

In our case the phase 5 describes an ordinary Josephson junction and the 
phase 9 an axion that passes through this Josephson junction. A very simple 
coupling scheme is given by 

5 + u 2 sin5 = j x (9 — 5) 



RC 



m 2 c 2 



9 + 3H9 + ~^sin9 = lx {5-9). (16) 



This corresponds to the case L —> oo in eq. (TT5j) . 

If the axion mass is at resonance with the Josephson plasma frequency, 
m a c 2 = hw, then synchronization effects of the phases 5 and 9 will oc- 
cur if 7a- is not too small, just as they occur for coupled Josephson qubits 
[j~9"l [20] . If the axion couples to fluxes similar as in eq. f|T5|) . then this is 
described by a small additional self-interaction potential given by V(a) = 
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— ^(|a 2 + 27rM 2 / a a). Quantum mechanically, one could even speculate on 
the formation of entangled states between axions and Josephson qubits. 

Given the quantitative agreement between the parameters of axion physics 
and Josephson junction physics outlined in section 3, one might hope that 
for axions the coupling ^ x is again of similar order of magnitude as in current 
nanotechnological experiments, provided the Josephson plasma frequency is 
close to the axion mass. This can be experimentally tested. 

There are no astronomical constraints on the size of ^ x since almost all of 
the matter in the universe is not in the form of Josephson junctions. Axionic 
dark matter may look completely 'dark' in the universe as a whole but not 
at all 'dark' in special superconducting devices designed by mankind, due to 
SQUID-like interactions. The only way to constrain j x is to scan a range of 
plasma frequencies and look for the possible occurence or non-occurence of 
universal resonance effects, produced by axions of the dark matter halo that 
hit terrestrial Josephson junction experiments. The intensity of this effect 
might display small yearly modulations, just similar as in the DAMA/LIBRA 
experiments pQ. What corresponds to tuning the cavity frequency in the 
experiments [8] would correspond to tuning the plasma frequency u in these 
new types of nanotechnological dark matter experiments. 

5 Axionic Josephson effect 

As an application of our theoretical treatment in the previous section, let us 
now discuss the analogue of the Josephson effect for axions, similar in spirit 
to what was experimentally observed for BEC in [25J. A Josephson junction 
biased with voltage V generates Josephson radiation with frequency 

ftwj = 2eV. (17) 

For such a biased junction the phase S grows linearly in time, i.e. 

2eV 

5(t) = 5(0) + —t. (18) 
h 

The relation between bias current / and applied voltage V is 

V = Ry/P - II « RI for / » I c . (19) 
Josephson oscillations set in if 

I > I c , (20) 
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i.e. the bias current / must be larger than the critical current J c of the junc- 
tion. In the mechanical analogue, the pendulum rotates with large kinetic 
energy. 

According to eq. ([2]) and (jSJ), the axion misalignment angle 9 will also 
start to increase linearly in time if it is being forced by very strong products 
of E and B fields. So from a formal mathematical point of view, an axionic 
Josephson effect is possible. The axionic Josephson frequency is given by 

huj = 2eV « 2eRI = ^^c 3 E ■ B, (21) 

where in the last step eq. (ITT]) was used. Condition ( 1201) translates to 

e 2 h 2 EB > flmlc = A 4 c (22) 

QCD-inspired models of axions fix A to be about 78 MeV Strong mag- 
netic fields in the laboratory correspond to about 10 Tesla, and strong electric 
fields to about 10 9 V/m. This gives E ■ B w 10 10 VT/m. One can easily check 
that under normal laboratory conditions one cannot produce stationary E 
and B fields of sufficient strength to satisfy eq. (T22|) . 

However, there is another interesting possibility how we can briefly induce 
axionic Josephson oscillations in the lab. This is based on the SQUID-like 
interaction mechanism discussed in section 4. Remember that the phase of 
a SQUID formed out of a Josephson junction and a passing axion should be 
gauge-invariant, provided the axion mass is at resonance with the plasma 
frequency. This led us to derive eq. ([TBI , meaning that the phase S of the 
Josephson junction and the phase 9 of the axion synchronize. For a bi- 
ased RSJ performing Josephson oscillations of frequency uj, synchronization 
means that 

5(t) = 5(0) + ujt (23) 

induces 

9(t) = 9(0) + ujt (24) 

for the axion. According to eq. (TTTT) . this means that the axion formally 
sees a huge product field E ■ B. The huge (virtual) magnetic field will make 
it immediately decay into two microwave photons. The microwave photons 
produced by axion decay have the frequency of the Josephson radiation and 
produce distortions in the I — V curve. They can be potentially measured 
in form of Shapiro steps (Shapiro steps are well-known step-like structures 
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in the I — V curves of irradiated Josephson junctions [2B]). Our theoretical 
idea thus opens up the possibility to develop new detectors for axionic dark 
mattei0. 

6 Conclusion 

Let us conclude. In this paper we have discussed the possibility that axions 
could interact with Josephson junctions by briefly forming SQUID-like inter- 
action states. Josephson junctions, SQUIDS, and spatially extended arrays 
of these superconducting devices can nowadays be built for a wide range of 
different parameters R, C, I c , and it is very easy to tune the bias current / 
to any value of interest. It is also very easy to adjust the plasma frequency 
of a Josephson junction to any value of interest. It is thus possible to build 
up a suitable resonant environment that could help to detect incoming dark 
matter axions. 

There is the prospect of developing new generations of detectors for dark 
matter axions that search for possible resonance effects and phase synchro- 
nization if the Josephson plasma frequency is close to the axion mass. In 
this way the size of the coupling 7 X between axions and a given Josephson 
environment can be experimentally constrained. As shown in section 3, the 
relevant dark matter mass parameter range is accessible by modern tech- 
nological developments in nanotechnology. An obvious advantage of these 
types of experiments is that the formal existence of extremely large prod- 
ucts of electric and magnetic field strengths EB acting on the axion can 
be simulated by a very simple experimental setup, an easily tunable bias 
current J, assuming that some axions hitting the Josephson junction will 
synchronize their phase due to a SQUID-like interaction. This effect may be 
systematically tested in future experiments. 
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